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Abstract 

Symbolic methods of umbral nature are playing an important and increasing role in the theory 
of special functions and in related fields like combinatorics. We discuss an application of these 
methods to the theory of lacunary generating functions of Laguerre polynomials for which we give 
a number of new close form expressions. We present furthermore the different possibilities offered 
by the method we have developed, with particular emphasis on their link with new family of special 
functions and with previous formulations, associated with the theory of quasi monomials. 
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I. INTRODUCTION 



It has been shown that symbohc methods of umbral nature provide powerful tools to deal 
with the properties of special polynomials and functions Q. These techniques greatly 
simplify the problems underlying such studies and allow to reduce the derivation of the 
relevant properties to straightforward algebraic manipulations. 

The key element of the game is the introduction of a symbolic operator used as an ordinary 
algebraic quantity. Within such a framework the ordinary Laguerre polynomials js, 4|^ 

Ln{x, y) = n! V / ^ — TT 
[rl)^[n — r)l 

are expressed as 

L„(x, y) = {y- cxYifo, (1.1) 

where c is an operator satisfying the property (fc^ = c"~^^, and acting on the "vacuum" ipo 
in such a way that 

1 

c (po 



T{l + a) 

The properties of Laguerre polynomials could have accordingly been reduced to those of a 
Newton binomial, containing an operator treated, in all the manipulations, as an ordinary 
algebraic quantity. 

We also introduce the family of associated Laguerre polynomials 

Ll^^\x,y) = ^-^^^^A^:\x,y) (1.2) 



with 

Al:\x, y) = nv - cxTvo = n\J2 u ^ ^ V ^^'^^ 

r\[n — r)!i (1 + a + r) 

Using the umbral method, the double- and triple-lacunary Laguerre polynomials can be 
expressed via aI^\x, y) as follows 

L2n{x,y) = {y- ax)2>o = {y- cxT \{y - cx)>o] = (-^) V^AW(a^,i/) (1-4) 

n \ / 



r=0 



^ Strictly speaking Laguerre polynomials are one- variable polynomials and are defined as L,i(.t) = Ln{x, 1). 
It is also easily checked that Ln{x, y) = y"Ln . We will use two-variable forms for future convenience. 
In final results we shall sometime use again the one-variable form. This comment applies to two-variable 
Hcrmitc polynomials, see Eq. (|2.4p . too. 
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and 



\3n. 



^0 = E t) it) (-^r'y'''-'-'^':^'\^^ y)- (1-5) 



L3n{x,y) = {y- cx) 

k,r=0 

A further family of polynomials, introduced by the same means as before, is provided 
by what we will call, for reasons which will be clear in the following, Laguerre- Wright 
polynomials, namely 

Ai" ''^) (x, y) = ny - c^xT = n\Y^ ^ (1.6) 

r\[n — r)!l [pr + 1 + a) 

From Eq. (O]) Al^'^\x,y) = Al?\x,y) follows. 

An interesting example for illustrating the flexibility and the usefulness of the symbolic 
methods is derivation of the generating functions of special polynomials and of special func- 
tions as well. We will consider in the following two types of generating functions, namely 
the exponential and the ordinary ones. In the first case, i.e. the exponential generating 
function of polynomials in Eq. (11. 6p leads to the following form 

EG^'^^^\x,y\t) = ^i^A^«(a;,2/) = ^-c"(2/-c^a;)" 

n=0 ■ n=0 

= e^*c"e-^%o = e^*PF('^'°+^)(-te), (1.7) 

where 

W/(/3.")(a;) = c"-i exp{c^x)^o. (1.8) 
In the case of ordinary generating function, Eq. (11. 6p gives 

oo oo 

oG^^'^\x,y\t) = 5^rA("'^)(x,y) = 5^rc-(2/-c'^x)" 

n=0 n=0 

1 c" 1 ^ f -tx \ 

tn:^fo = — 7-Ep,a+i - — - , (1.9) 



l-tyl + ^c/^"'^ l-ty ^^-^^\l-ty 
where 

Ep,^{x) = T^^^^o. (1.10) 

1 — cPx 

In Eqs. (II. 8p and (ll.lOp we have formally defined the Bessel- Wright and Mittag-Lefiler 
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functions, respectively js-T], which in expUcit form read 

oo oo 



X 



r=0 



r=0 



-^0 



X 



r\T{(3r + a) 



X 



T{(3r + a) 



r=0 r=0 

The generating function of the associated Laguerre polynomials are accordingly derived 
as the special case of previous demonstrations: 

oo oo n 

oGS"'^)(x,y|t) = ^r4")(x,|/) and ^gS"'^)(x, = ^ y), 



n=0 



n=0 



oG^r''\^,y\t) 



-l-a 



ii-tyy+^ 



n=0 



exp 



(1 - tyy^- ^ 



(1 - ty) 



1+a 



(1.11) 



which for lI^\x) = L\^'{x^l) gives the well-known formula (8.975.1) on p. 1002 of 



Furthermore, repeating similar calculation for exponential generating function for integer 
a = m, m G N, we have 

{n + m)\ 



n=0 



■[t{y-cx)r 



r=0 
r=0 



.s=0 ^ ^ 
o— n ^ / 



m \ m! 



m\m\ 
V 1 t\ 



t)W(^'™+'+^)(-ta;) 



(1.12) 



where we used the binomial formula and the well-known relation 



(i-^)"" = E 



n=0 



r(n + a) 
^ r(a) ^ 



\z\ < 1. 
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II. LACUNARY GENERATING FUNCTIONS OF LAGUERRE POLYNOMIALS 



The double-lacunary exponential generating function of ordinary Laguerre polynomials 



nl 



(2.1) 



n=0 



is apparently not known. A related expansion given in literature does not appear to be 
correct, see Eq. (5.11.2.10) p. 704 of {q]. In the following we obtain the explicit form of Eq. 
(12. ip in terms of known functions. 

According to our procedure, we rewrite Eq. (12. ip using Eq. (I1.4p as 



^-^ nl 



cx 



,2n 



n=0 



= e^'*e^'"'*-22^^"Vo- (2.2) 



We can provide a definite meaning for the previous expression in terms of known special 
functions by recalling the following expansion 



bz'^+az 



'yTl 



(2.3) 



n=0 



with 



where [n 
mials 



[n/2] 



Hi'\a,b)=n\Y^ 



r=0 



(n — 2r)!r! 



(2) 

the floor function. The polynomials Hn (o^) are two- variable Hermite polyno- 
10|, also defined through the operational rule |ll| 



They reduce to the ordinary Hermite polynomials Hn{z) through the relation 



i/f (x,y) = (-^y^)"i7„ 



IX 



(2.4) 



The above family of polynomials provides a basis for the definition of the so-called if-based 
functions [ll; for example the if-based cylindrical Bessel functions are defined as 

(2) 



r=0 



They have been obtained by replacing x""'"^'' in the relevant series expansion by H^^2ri^i v)- 
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According to Eq. (12. 3p . we obtain 



r=0 



2^ /,,n9 HriiyxVt) 



r=0 



(2.5) 



where 



J? 



C*n. (X. 



H^'\-x,y) 



n ^'r(i 

r=0 ^ 



!r(l + a + r) 



X] 



is the if-based version of the Bessel-Tricomi function 

oo 

Ca{x) = ^ 

Thus we get in conclusion 

.(1,1) 



r=0 



r!r(l + a + r) 



sG'^'''^(x,y|t) = e^ *^,Co(-2x2/t,x2t). 



(2.6) 



In terms of standard Laguerre polynomials Ln{x) = Ln{x, 1) Eq. (12. ip can thus be rewrit- 
ten as 



oo 



oo 



—L2n{x) = e* V 



(2.7) 



n=0 ■ ■ r=0 

which constitutes one of key results of the present investigation. The validity of identity ()2.7p 
as well as of further results of this type for other generating functions (see Eqs. (I1.12p . and 
(12. Sp . (I2.14p . (13. ip . (13. Sp below) can be independently proven by the substitution t —t"^, 
followed by coefficient extraction and the use od Eqs. (II. lip and (12. 3p . see |l2 |. 

Using the above technique, see Eqs. (I2.ip -( l23|) . we derive the exponential double-lacunary 
generating function for associated Laguerre polynomials Ln"''(x), for a = m 
and establish that 



1,2, 



n=0 



r=0 



r!(r + 3m)! 
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which for standard associated Laguerre polynomials L^\x) = L^\x,l) and using Eq. 
(12 ■4p has the form 

P2m{r;x, l,t) 



J. '. „ 



r!(r + 3m)! 



[tx 



VtyHr{iVt). 



(2.8) 



The polynomials p2m{f", x,y,t) are of degree 2m in the variable r, with the coefficients 
depending on x, y and t. For m = 1, 2, p2m{r; x, y, t) have the explicit form 

P2{r; X, y, t) = {1 + 2ty^)r^ + (5 - 4:xyt + 10ty'^)r + (6 + 12ty^ - Uxyt + 2tx^), 



Pi{r; X, y, t) = (2 + lOty^ + AtY^ + [36 + (ISOy^ - 20yx)t + {72y^ - 16y^x)t'^y 



+ [238 + (lOx^ + 1190?/2 - 300yx)t + {-240y^x + 24y'^x^ + 476?/^)t>^ 



+ [684 + (llOx^ - USOyx + 3420y^)t + {-llSAy^x + 264?/V - 16yx-^ + 1368?/^)r]r 



+ [720 + (-2400?/x + 3600^/2 + 300^^)^ 



+ (-1920y^x - 9Qyx^ + 720yV + 1440y^ + 4x^^)1^] 



With a moderate effort the polynomials P2mij] y, t) for m > 2 can be also obtained. 
The method can be extended to the more general cases like 



EGi''^'\x,y\t) = J2'-Aif\x,y), 



n=Q 



which, according to our procedure, can be written as 



r=0 



EG^t^\x,y\t) = c^y-{y-c^xy^ipo = ey"'c''e-^''^y'+'"''^"'ipo 



r=0 



^0 



^/'HW^f'^+^\-2xyt,xH) 



(2.9) 



with 



H 



W^^'^\x,y) = f2^p^. 

^-^ r!i [ar + p) 



r=0 



being the if-based version of the Bessel- Wright function. 



We give here without demonstration the following two formulas of similar type 



n=0 r=0 



oo (r) 



(1/2). 



tx 



2(1 -t) J 



and 



n=0 r=0 ^ 



r! — X 



(r-s)!(r + 2s)! " V3 



(2.10) 



(2.11) 



In this section we have shown that the tools of employing the symbolic method to de- 
rive lacunary generating functions for Laguerre polynomials can easily be applied with a 
minimum of computation effort. 

We now employ the symbolic method to derive old and new generating function for the 
associated Laguerre polynomials. We start with simple observation, namely 



c e </5o = > ^ -r*^ V^o ^ 

n=0 



(2.12) 



r(l + a) 

n=u 

Thereafter, we calculate Q[°'\x,y\t) which after applying Eqs. (11. 2p . (I1.3p and f l2.12p has the 
following form 

^ t I y — cx 

n=0 \_n=0 ' ^ 

= e-*"r(l + a)c" ey'^'" </^o = (1 + Z/t)"e-*", (2.13) 



g[''\x,y\t) = X^r4"-")(x,y) = c°r(l + a) 



Formula fl2.13p for the ordinary associated Laguerre polynomials, i. e. for y = 1, is equal to 
(5.11.4.8) on p. 706 of {q]. In the case of generating function 

oo 

gi''\x,y\t) = J2eLt'''\x,y) 

Eqs. ([L21) and ([L3l) yield 



n=0 



e^°'(rc,j,|i) = c"r(l + a) 



.n=0 



f y — CX 



(2n)! 



2n 



C°T(1 + a) cosh ( — X\ft I <y9o 



r(l + a) 



cosh(a;A/t) cosh ( -^^C^^ — sinh(xA/t) sinh ^ 



[1 — ty"^) 2 |cosh(A/tx) cos(T) — i sinh(\/tx) sin(T) | 
1 -t?/2)tcosh(v^a;-«T), 



(2.14) 
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where 



T{t, y) = a arcsin 



III. FURTHER DEVELOPMENTS 



The method we have developed so far can be extended to obtain shghtly more comphcated 
expressions hke 



= —,L2n+i{x,y) and = —L^n+i{x,y), 

n=0 n=0 

for / = 1, 2, . . .. Note that corresponding formulas for the Hermite polynomials were obtained 
in Applying Eqs. fll.ip . (12.21) . and (12. 3p for the generating function of double-lacunary 
Laguerre polynomials, we have 

I 



,G^}i'\x,y\t) = ey'\y-cxy 



r=0 



Y,-Hi'\-2yxt,xH) 



(po = e 



s=0 



E W-\-xrHCs{2xyt,xH), 



which, when written in terms of standard associated Laguerre and Hermite polynomials, is 
equal to 

I 



y^—L2n+i{x) = e' 
^ — ' n\ ^ — ' r 

n=0 r=0 



oo 



. v-^ [ixy/tY , . r-. 



s J [r + s) 



r=0 



r\[L + r)\ 



Using the analogous procedure we obtain also EG^i'f'i'^i y\^) form 



EGtfixMt) = ey'\y-cx)' 



J2^Hi:'\-^^y't,3x'yt,~xH) 



r=0 



y [-x)h 



Ci^\3xyh,3x^yt,xh), 



where 



HCi^\x,y,z) 



°° Hi^\-x,y,-z) 

r=0 



r\(r + s)! 



is a third order if-based Tricomi function, with 



r=0 



(n — 3r)!r! 



(3.1) 



(3.2) 



being a third order three variable Hermite polynomial, with generating function [l^ 

n=0 

Eq. fl3.2p expressed via standard Laguerre and Hermite polynomials is given as 



^-^ n\ ^-^ in 



n=0 



71=0 



n + l)\ 



[n/3] 



-J7Z —, ^n-3r I i 



r=0 



r\(n — 3r)! 



(3.3) 



In the case of the associated Laguerre polynomials for / = 0, L^^J^{x,l) = L^J{x), the 
generating function is given as 



E V3n(^) = e*V^^^^iJ,,(-3tx,3te^-tx 



q3{r;x,t) 



n=0 



r=0 



E 

r=0 



r!(r + 4)! 



q3{r;x,t) 
r!(r + 4)! 



[n/3] 

E 

r=0 



{-tx^YiixV^)''-^'' / V3t 

, , . I J. J. St" I 

r!(n — 3r )! 



(3.4) 



with 



g3(r; t) = (1 + 3t)r3 + (9 + 27t - <dtx)r'^ + (26 + 78t - 63tx + 9^^)^ 

+ (24 + 72t - lOStx + 36^2 - 3te^), 

which is a third order polynomial in r, with the coefficients depending on x and t. The 
extension to the case 



n=0 



can be straightforwardly accomplished within the present framework and writes 



E V^n-«i,«2,...,(-l)'"«. 



,r=0 



m 



s=0 



ai, ^2, . . . , amj, 



V 



x'Py'^-n, p=l,2,...m, 



(3.5) 



where hC^^ (ai, 02, ... , a^) is an m-th order based Tricomi function defining as basis the 
m-th variable Hermite polynomials H^\—ai, a2, ■ ■ ■ , {—ly^am) specified by the generating 
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function 



—H^^ {xi,X2, . . . , Xm) = exp I ^ Xst"" 

n=0 ^' \s=l 

Let us finally consider the bilateral generating function 



(3.6) 



^ -j- 

—Lnix,y)Lniz, 
„ n! 



u\ 



n=0 



t 



n=0 



</2l,0V^2,0 



(3.7) 



r(l + a 



■, g = l,2. 



Accordingly we obtain 

E— L„(x, y)Ln{z, u) = e*"^ exp {-cixut - C2yzt + ciC2X2;t) v5i,oV52,o 

= e^^^Cofi{—xut, —yzt\xzt), 



n=0 



(3.8) 



where 



//(^o,o(a;,l/|r) = ^ 



r\slk\{r + k)\{s + k)\ 



r,s,k=0 

Closely related considerations devoted to combinatorics of Laguerre polynomials are devel- 
oped in I15I ]. 

Using the operational method we can also obtain the other interesting identities of the 
generating function of the lacunary associated Laguerre polynomials. These formula are 
listed below: 



^ (l + -) 

V ^ 2 A 



^iJ (^) 



00 Y (^+0) 



(!) 



^ (l + -) 

K"- ^ 2Jr 



tx 



2(1 -t) J 



00 n\ 



n=0 



(1 + m)n 



x^/t\ ^ ( tx . ^ 
exp I -Y~t 



X 



Vt 



1-t 



(3.9) 



(3.10) 



where m = 0, 1, 2, . . . and Im{z) is the modified Bessel function 

^ (i + f)„(l + f), ' ' + - 



E 

.5=0 



,.„(i + f ).(! + §) 

(f) 

V-s)! r(2s + a + r + l) 



tx 



9(1 -t) 



(3.11) 



Eq. ( KTO^ corrects the formula (5.10.2.10), p. 704 of |9|. 
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Before concluding this paper it is worth presenting some further comments to reconci 



the present results with previous approaches based on the monomiality principle 



16 



-ll8|. 



Within that context the so called Laguerre derivative lD^ has been introduced, so that 

iD^Lnix, y) = nL„_i(x, y). (3.12) 

The Laguerre derivative in differential terms has been defined as = —dxxdx- It is evident 
that, according to the umbral technique proposed in this paper, the Laguerre derivative can 
alternatively be introduced as 

lDx = -c-'dx. (3.13) 

On the basis of such a definition we find 

g?/L^.g-x _ e-y£-'d.^_^^ = ^ = _J: 

1 + cx 1 + c{x — yc~^) 1 — y 1 -\- 



thus recovering the well-known result 18 1 



1 - y 

As a further example of application we consider the action of the exponential containing the 
Laguerre derivative acting on the Tricomi function, namely 

e^^^-Co(x) = e-^"^"'^-e-'^>o = e^Co(x), (3.15) 



in agreement with the fact that Co(x) is an eigenfunction of the Laguerre derivative 16 1. 
The formal procedure we have envisaged can even be pushed further. 
The use of the following identity 

K^9-f{x) = f{Kx). (3.16) 

and the definition of the 0-th order cylindrical Bessel function as a pseudo- Gaussian [l|, 
namely 

Mx) = e-<t)Vo (3.17) 
allow the derivation of the following identity 

c-^^'^^Jo(x) = e-(t)', (3.18) 



which, in terms of integral transforms, can be interpreted as a kind of Borel transform [i 

roo 

c-5^^-Jo(x) = y e~'Jo(^Vsxjds. (3.19) 
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Finally let us note that having expressed the cylindrical Bessel functions in terms of Gaus- 
sian, it is also possible to "reduce" a Gaussian to a Lorentzian, according to the identity 

1 + cx^ 

According to this last identity we can write the relevant integral as 

J e~^^d^ = c~2 arctan (^Vdx^ipo. (3-21) 

Here we give for completness the list of new close form expressions obtained in the present 
investigation. These are equations: ([L7D, (US]), (HHD, dM]), (EH), (jM]), ^AO\f . (I2ll|) . 
(EUD, dSH), (E3D, (IS3D, (IS3D, (E3D, f lXTU]) . and fIXTT]) . The methods we have illustrated in 
this paper appear fairly flexible and amenable for further implementations, as will be shown 
in a future investigation. 

We thank A. Bostan for important discussions. This work has been supported by Agence 
Nationale de la Recherche (Paris, France) under Program PHYSCOMB No. ANR-08-BLAN- 
0243-2. 
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